A quantum topological phase transition at the microscopic level 
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We study a quantum phase transition between a phase which is topologically ordered and one which is not. 
We focus on a spin model, an extension of the toric code, for which we obtain the exact ground state for all 
values of the coupling constant that takes the system across the phase transition. We compute the entanglement 
and the topological entropy of the system as a function of this coupling constant, and show that the topological 
entropy remains constant all the way up to the critical point, and jumps to zero beyond it. Despite the jump in 
the topological entropy, the transition is second order as detected via local observables. 



I. INTRODUCTION 

Some strongly correlated quantum many body systems dis- 
play a type of order which cannot be characterized by any 
local order parameter Instead, such order is topological in na- 
tural, with the fractional quantum Hall systems being the pri- 
mary example so far: they are liquid states that exhibit exotic 
properties such as a ground state (GS) degeneracy that can- 
not be lifted by any local perturbations^'^* and fractionalized 
degrees of freedom"^. Topologically ordered states are also 
interesting in that their robustness against local perturbations 
might be of use for decoherence-free quantum computation^. 

An example of an exactly solvable lattice spin model that is 
topologically ordered was presented by Kitaev in Ref. 5, and 
the system was argued to be robust against small perturba- 
tions that tend to order the system a la Landau-Ginzburg and 
take it away from its topological phase. The departure from 
the topologically ordered phase should occur through a quan- 
tum phase transition. Such quantum phase transition, how- 
ever, cannot be entirely captured by ordinary methods based 
on local Landau-Ginzburg order parameters, and new methods 
need to be devised in order to investigate the fate of topolog- 
ical order across the phase transition. These novel methods 
must be based on the fundamental properties of topologically 
ordered phases, such as the GS degeneracy in presence of a 
gap, and the presence of a non-vanishing topological entropy. 

Recent efforts to understand quantum phase transitions in 
topologically ordered states include a mean-field approach for 
these exotic states*, and analytical and numerical studies^'^''' 
of the Kitaev model in the presence of a field. The numer- 
ical analysis presented in Refs. [^Hl leads to the conclusion 
that topological order survives unchanged up to the second 
order phase transition at — 0.32847(6) (in the notation of 
Eq. (fTSll). while the system is no longer topologically ordered 
for (3 > j3^. (Here (3 stands for the coupling constant that 
drives the T — Q quantum phase transition - the notation will 
become apparent shortly, and is chosen because of a close re- 
lation to a classical model.) 

In this paper, we investigate analytically a quantum phase 
transition out of a topological phase. We show that the re- 
cently defined topological entropy'" " works well as an "order 
parameter" across the transition. We study the transition using 
a model - see Eq. (|5]l - that is shown to behave much like the 
Kitaev model in a magnetic field for small values of the field. 



The advantage of this model is that the ground state can be 
obtained exactly, from which we can then compute the topo- 
logical entropy explicitly, and show that it remains constant in 
the topologically ordered phase (/3 < ~ 0.4406868), drop- 
ping abruptly to zero in the non-topologically-ordered phase 
(/3 > (3^, despite the continuous (second order) character of 
the transition. 

We find that in this model, even though one cannot iden- 
tify a local order parameter that vanishes in one phase and not 
in the other, one can show that the (continuous) local magne- 
tization has a singularity at the critical point. In the model, 
we show that the magnetization equals the energy £'ising(/3) 
of a 2D classical Ising model with N spins evaluated at an 
inverse (classical) temperature equal to the value of the cou- 
pling constant (3 that drives the system through the T = 
phase transition: 



(1) 



From this relation, it becomes evident that the magnetization 
m{(3), although continuous and non-vanishing across the tran- 
sition at 13c (much as the energy of the classical Ising model 
across the classical transition), has a singularity in its first 
derivative, since 



dm 
'8(3 



1 dEising 

N d(3 



(2) 



and the Ising model heat capacity Cising diverges logarithmi- 
cally at j3c- Hence, although there is no local order parameter 
that can detect either the topological or the non- topological 
phase in this system, one can expose the topological quantum 
phase transition through the singularity in the derivative of a 
local quantity. 

This is contrasted, for example, with the case discussed in 
Sec, nil B1 where a similar topological transition is accompa- 
nied by a simultaneous Z2 symmetry breaking phase transi- 
tion. In that case, the same transition is captured both by 
the non-local topological entropy, and by a local (Landau- 
Ginzburg) order parameter. 



II. THE MODEL 

The model that we consider is a deformation of the Kitaev 
model^, and it is defined on a square lattice with spin- 1/2 de- 
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grees of freedom living on the bonds, as shown in Fig. [T] The 
pure Kitaev model is written in terms of star and plaquette 
operators (see Fig.[TJ. Star operators are defined as 



n ^-11' 

iGstar(s) iSs 



(3) 



where i labels the four spins on the bonds departing from some 
vertex s of the square lattice. Plaquette operators are defined 
as 



B„ 



n ' 

2Gplaquette(p) 



(4) 



where i labels the four spins on the bonds around some pla- 
quette p of the square lattice. 
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FIG. 1: (Color online) - Examples of star and plaquette operators, 
centered at a lattice site s (blue open circle) and at a dual lattice site 
p (red open square), respectively. The solid black dots represent the 
spin-1/2 degrees of freedom living on the bonds of the lattice, and 
the dashed lines connect the spins involved in the definition of each 
of the above operators. 

The Hamiltonian we consider in this paper is 



H 



= ^^Etaev + Ai^e-^S..,*^^ (5) 

s 

where Ao,i > and is a parameter that we use to tune the 
system across a topological quantum phase transition. Notice 
that for f3 — the Hamiltonian Eq. (|5]) is simply the Kitaev 
Hamiltonian -ff kitaev Ref 5, up to a trivial overall constant 
shift of the energy. 

The exact ground state wavefunction of this Hamiltonian 
can be obtained by deconstructing H into two pieces, H = 
XqHq + XiHi, as follows. 

Take G to be the (Abelian) group of all spin flip operations 
obtained as products of star type operators. Notice that g'^ = 1 



for any element g of the group G. By acting with elements of 
G on a given reference configuration \a\) one generates 
a manifold of states, which however does not encompass the 
whole basis. For example, the action of a star operator can- 
not change the sign of the product of ct^'s around any square 
plaquette in the lattice^'^^ (see Fig. [1]). Therefore, there is a 
non-trivial (and non-unique) minimal set { | ^) } of reference 
configurations that generates the full (T^-basis under the action 
of the group G. (In particular, one of the elements in this set is 
the reference configuration |0) that is fully magnetized in the 
z-direction, say erf = 1, Vi.) 

Consider then the family of Hamiltonians 



i?i(/3) = E 



(6) 



for some real-valued parameter (3. The ground state of any 
such Hamiltonian can be obtained exactly and it can be written 
in the form 

\GS,) = Y.^o.Y. 7^ (7) 



see 

where a labels the different block-diagonal sectors corre- 
sponding to the states in the minimal set {IvP^)}; af{g,a) 
is the z-component of the spin at site i in state g I^'q,); and 
the coefficients can be chosen at will, subject to the nor- 
malization condition IV-'qP = 1- Although the choice of 
minimal set is non-unique, one can show that Eq. (|8]l is in- 
dependent of such choice, modulo an irrelevant permutation 
of the a indices. Within each block-diagonal sector, the GS 
of Eq. (|6j is unique. Instead of proving directly that ([8]l is 
the GS of (|6|l, it is more convenient to notice that the family 
of Hamiltonians in Eq. (|6]l is a particular choice of Stochastic 
Matrix Form decompositions of quantum Hamiltonians that 
exhibit precisely Eq. ^ as their GS.'^ One can verify this by 
showing that each of the operators 



Qs 



(9) 



between square brackets in Eq. (|6]l annihilates the inner sum 
in \GSi), independently of the index a. That the GS energy 
is zero follows because 



and 



, cosh /3 E a-," 



0, 



(10) 



(11) 



from which it can be shown that the expectation value of Qs 
with respect to any state is always greater than or equal to 
zero. 
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Let us consider now the remaining part of the Hamiltonian, 



One can use the decomposition 



^0= E n 

P 'iGplaquette(p) 



p i&p 



(12) 



Recall that any star operator A^, and therefore any element of 
the group G, preserves the product Yliep '^^ every plaquette 
of the lattice. The GS wavefunction of Hamiltonian (fT2] i can 
then be written as 



\GS, 



a g^G 



(13) 



for any choice of the coefficients 0^ „ (J^J J2geG I<?^s,qP = 
1). Here the primed sum over a is restricted to the (four) 
block-diagonal sectors that satisfy Yliep ^\ ~ ^'-'^ P^^" 
quettes p in the lattice, and it must be carried out separately 
because no operation in G allows to change sector^ 

As a result, any linear combination with positive weights 
Aq and A^, 



= ^KUaev + AiE^"^^' 



(14) 



and therefore our Hamiltonian in Eq. (|5]l, has the GS given by 
|G5)=E^"E 7^ 3l*">- (15) 

a gee 



Notice that one of the topological sectors that satisfy 
riiGp '^^ ^ "^P '^^^ containing the fully magnetized 
configuration in the z-direction (|0)). 
In particular for ^ 1 



e-/3 5:.ee*? ~ (16) 
AiE^"'^^'"''*' - const - 2/3Ai ^ CT^ (17) 

s i 

Therefore, in the limit of small f3 (in absolute value) the 
Hamiltonian in Eqs. ( I5I14| | is equivalent to the Kitaev model 
in presence of a magnetic field proportional to /3A]^, 



H 



-^0 E n '^^^ - ^1 E n -^^^ - 2/3A, E 

p i£p 



(18) 



ies 



For larger values of /3, the many-body terms in Eq. ( fT9] l are 
no longer negligible and the equivalence is lost, although the 
form of the GS (fTsT i suggests that the system gets deeper and 
deeper into the magnetized phase - as one would expect upon 
increasing the strength of the magnetic field in the Kitaev 
model. As we discuss in Section 1111 Al our model undergoes 
a second-order phase transition at = (1/2) ln(\/2 + 1) ~ 
0.4406868, where it displays a dimensionality reduction that 
places the transition in a different universality class than the 
one studied in Refs. J,i8i^ 



e-/3E.Ee*t ^ [| [cosh(/3) -CT,"sinh(/3) 

= cosh''(/3) 

- cosh3(/3)sinh(/3)E'^i 

+ cosh2(/5)sinh2(/3) ^ ct^^ct^^ 

- cosh(/3)sinh3(/3) 

+ sinh4(/3)[]a^ (19) 

ies 

to estimate the limit of validity of Eq. ( fTsT i to be given by the 
condition 



cosh3(/3) sinh(/3) 




cosh(/3) 


cosh2(/3) sinh2(/3) 




sinh(/3) 



> 2. 



(20) 



This corresponds to a ratio between the coupling to the mag- 
netic field and the coupling to the cooperative transverse field 

(E. n»es ^Y) 



2Ai cosh^(/3) sinh(/3) 



< 



16 



1.78. 



(21) 



The detailed numerical analysis presented in Refs.|8||9|lead the 
authors to conclude that topological order survives up to the 
second order phase transition at finite /J^ (in the notation of 
Eq. (fTSl)), while the system is no longer topologically ordered 
for /3 > /3^. In the following, we investigate this phase and 
the relative phase transition using the exact ground state of 
our model (|5]l to compute the topological entropy'" " across 
the transition. Using an exact derivation from the microscopic 
degrees of freedom, we show that the topological entropy is 
able to detect a transition from a topologically ordered phase 
(fi < j3^) to a non-topologically-ordered phase (/3 > (3^. In- 
deed, it remains constant at its known /3 ^ value up to the 
transition and drops abruptly to zero afterwards, despite the 
continuous character of the transition. 



III. THE TOPOLOGICAL ENTROPY OF 
FACTORIZABLE (LOCAL) WAVEFUNCTIONS 



Using the definition in Refs. fTolfTll the topological entropy 
can be obtained as a linear combination of Von Neumann en- 
tanglement entropies of different bipartitions of the sys- 
tem into subsystems A and B: 



S'4 = -Tr [pa log2 Pa] -^vn- 



(22) 



where pj^ = Tr^ [p) is the reduced density matrix obtained 
from the full density matrix p by tracing out the degrees of 
freedom of subsystem B, and the last equality holds when- 
ever the full density matrix p is a pure-state density matrix. 
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The different bipartitions are aimed at removing all the exten- 
sive (boundary) contributions to uncover the sole topological 
contribution. A particular choice of the four bipartitions'" is 
illustrated in Fig.|2] and the topological entropy is then defined 



boundai-y 2 



boundary 2 




(2) 



(3) 



(4) 



FIG. 2: Illustration of the four bipartitions used to compute the topo- 
logical entropy in Ref.fTol. 



as: 



r,H — >oo \ 



(lA) 



(2A) 
VN 



c(3A) _ MiA) 
>JVN '-'VN 



(23) 



In order to compute the topological entropy as a func- 
tion of the parameter /3, let us first notice that the Hamil- 
tonian (fl4] i constructed above, with GS given by (flST l, be- 
longs to a class of Hamiltonians whose GS wavefunctions 



geG ' 



-PE/2 



g\0) have non-negative, fac- 
torizable amplitudes, i.e.. Eg = E^^ +E^^, with g = 9a®9b 
for all bipartitions [A, B). For this type of Hamiltonians, one 
can compute the entanglement entropy as follows. 

Consider a given bipartition {A, B) of the system. The re- 
duced density matrix = Tr^ (p), obtained by tracing over 
all degrees of freedom in B, is given by^ 



Pa 



- E 



X {0b\9b9b9b\0b) 9a\0a){0a\9a9a 

e ' "Be 

X 9a\0a){0a\9a9a 



1 \- 

z 



(24) 



where 9^ 9a® 9b^ |0> = \^a) ® |0b> and c G (G^ c 
G) is the subgroup of transformations acting solely on A (B) 
and leaving B (A) invariant: 

Ga = {.9 e G I 5b = Is} 
Gg = {geG\gA = 1^} 

Notice that we used the group property to rewrite a generic 
element of G as gg, 3! g € G, as well as the additive property 
oiEg. 

We can then compute the trace of the n-th power of the 
reduced density matrix Tr [(/O^i)"] and use the identity 

d 

- lim T-Tr [(p^)"] = -Tr [p^ In p^] (25) 

n-^l (jn 



to obtain the Von Neumann entropy = — Tr [p^ log2 



X e ' - ^'-^ <'i.A"^.A' 

X {0a\9ua9'i,a92.a\^a){^a\92.a9'2.a93,a\^a) ■ ■ ■ 
X •■•(0aI5„,a5'i,a5i,a|0a> 

g'i^---,g'n'^G^ 

X (0aI3i,a5'i,a52,a|0a>(0aI52,a52,a53,a|0a> • ■ • 

X ■■■{0A\9n,A9'n,A9l,A\^A)^ 

(26) 

where we used the fact that the inner products in Eq. (|26] | im- 
pose 



9i+l,A — 9i,A9i,Aj 



(27) 



for i = 1, . . . , n, with the identification n + 1 = 1, and there- 
fore 

The condition in Eq. ( |27] i can be satisfied if and only if 

9^ = i9^+l.A9^^A) ® Ifi ^ Ga- (29) 

Thus, the summation over all g' of the inner products in 
Eq. ( |26] | yields a constraint over the allowed values of 
ffi, . . . ,5„ e G: 

«) 1b e G^, = 1, . . . , n (n + 1 = 1) 

{9^ A9j a) ^B ^ G^, Vi,j = l,...,n 
9i9j^GaGb, Vi,j = 1, . . . ,n. 

In particular, the last line can be recast as 

9i = hi 9i fc, 3! hi £ G^, fc, e G^, V i = 2, . . . , n. 

(31) 

The physical meaning of these conditions will become clear in 
the next section, for the specific case of the system considered 
in this paper, although the form of Eq. ( [3T1 ) already suggests 
that they require all the elements g^ to agree at the boundary 
of the bipartition {A, B). 

We can finally use Eq. ( [3T| i to simplify Eq. ( |26] l 



(30) 



Mp1\ = jr.T.--'''^ E 



see 



h^,...,h^eGj^ 

A,Q 1 ' ' ' 7 ^ ^ ^ G TD 



z ^ 

geG 



eG.^keG, 



n-1 

I (32) 
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and obtain, via Eq. 



ff6G 



eG.,keG 



A' "-^^B 



l0g2 



Z 



= -(log2Z„)+log2Z = /3 (F„)-F , (33) 



where Fg is the partial free energy given by all the configu- 
rations that can be obtained from g via products of spin flip 
operators that act solely on subsystem A or subsystem B (i.e., 
having the same 'boundary' as g), and (. . .) denotes the en- 
semble average over g £ G with weight e^^^a . Notice that 
our result in Eq. (l33l l is the lattice equivalent of the Von Neu- 
mann entropy obtained by Fradkin and Moore in Ref. 14 for 
continuous systems. 

Alternatively, Eq. (|33l l can be interpreted as the entropy of 
mixing (or configurational entropy) of the allowed bipartition 
boundaries in G. This can be made more transparent by intro- 
ducing the quotient group Q = G/ (Gj^Gg), and by rewriting 
Eq. ([33ll as 



5, 



(A) 



EE 



X loe 



E^« l^Ss^' 



z 



(34) 



where we used the fact that the term in square brackets is in- 
dependent of h and k, and where we introduced the notation 



P„ 



E 



heG.,keG^ 



z 



(35) 



for the probability of boundary q to appear in G, for a given 
inverse temperature /3 and energy . 

In order to proceed further, let us focus for simplicity on the 
specific GS of our system (flST l. The generic case of a wave- 
function with factorizable amplitudes can be inferred with mi- 
nor modifications. 



A. The case of 1-body potentials 

All of the above results apply straightforwardly to the GS 
in Eq. ( fTSl l. Notice that (i) the topological entropy in each 
block-diagonal sector of the pure Kitaev model is the same^^ 
and (ii) it is reasonable to make the working assumption that 
the relevant sector for the transition to the fully magnetized 
state |0) is the one that contains this state, and that is therefore 
obtained upon applying the group G to |0). For the purpose 



of computing the topological entropy, one can thus replace 
Eq. dnil by 



|G5) = ^5:e'^2:.^K.)/2,|0). 



(36) 



and obtain 



geG 



X Ids 



E 



heG,,keGf 



aPT.,<^\{hgk) 



(37) 



where Z = X;geGe^^''^'^^ 



In order to simplify Eq. dJTl i with the purpose of comput- 
ing the topological entropy of the system ( |23] |. it is conve- 
nient to do the following change of variables. Recall that a 
generic configuration g|0) is uniquely specified by the set of 
star operators acting on the reference configuration |0), which 
we chose to be the ferromagnetic state with all the a spins 
pointing up, modulo the action of the product of all the star 
operators (which is equal to the identity). Thus, there is a 1- 
to-2 mapping between G = {.g} and the configuration space 
<d = {9} of an Ising model with degrees of freedom 6^ living 
on the sites s of the square lattice, where for example 9^ = —1 
(+1) means that the corresponding star operator is (not) acting 
in the associated g. Since each cr spin can be flipped only by 
its two neighboring 6 spins, then <t^ = 9^6^,, where i labels 
the bond between the two neighboring sites (s, s'), and 



E- 

geG 



eee 



(38) 



Notice that, using the above mapping, the GS wavefunction 
of our model, Eq. (fTsT l. can be rewritten as 

\GS) = ^ 9{e) |0), (39) 

eee 



where Z = J2eee ^ 



Thus, all equal-time cor- 



relation functions that can be expressed in terms of the 9^ 
variables are the same as those of a 2D classical Ising model 
with reduced nearest-neighbor coupling J /T = (3, implying 
that the critical point of the latter (3^ = (1/2) ln{^/2 + 1) ~ 
0.4406868 corresponds precisely to the critical point of our 
quantum system. Notice also that the magnetization in the 
original a spin language is indeed the nearest-neighbor spin- 
spin correlation (i.e., the energy) in the 9 spin language. 



-E 



eee 



eee 



^E-?(.9W) 



-E' 

{s,s') 



N 



-Rising (/?) 



(40) 
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Therefore, one concludes that the magnetization m{(3) is con- 
tinuous across the transition at f3c but there is a singularity in 
its first derivative 



dm _ J_ dEising 
dl3 ~ N dp 



(41) 



as the Ising model heat capacity Cising diverges logarithmi- 
cally at Pc- 

In the following, we will show how such continuous phase 
transition is accompanied by a sudden, discontinuous vanish- 
ing of the topological entropy of the system. 

The case of a configuration of the form hgk, with h e 
and k £ Gg requires a few additional steps. First of all, notice 
that the composition of any two elements g,g e G is repre- 
sented in the 9 spin language by the site-by-site product of 
the two configurations corresponding to g and g, respectively: 
es{9~9) = ^.(5) OM- In particular, e^{hgk) = 0,(/ifc) OM- 

Moreover, using similar arguments as in Ref. 15, the star 
operators of a bipartite system [A, B) can be divided into bulk 
star operators, i.e., those acting solely on subsystem A or sub- 
system B, and boundary star operators acting simultaneously 
on A and B spins. The boundary star operators can be fur- 
ther subdivided into different sets according to the different 
boundaries around each connected component of A and B (for 
a total of + — 1 boundaries, and to^ being the 
number of connected components of A and B, respectively). 

Let us define a collective operation as the product of all 
the star operators in one of these sets. That is, the product of 
all the stars around a connected boundary of the bipartition 
[A, B). Clearly, the number of such collective operations is 
given by the number of sets, to^ + rrig — 1. 

One can show that the subgroup Gj^Gg C G, to which 
the product hk belongs, can be generated by all the bulk star 
operators together with all but one of the collective operators 
(all but one, independently of which one is chosen to be left 
out, is requked because the product of all boundary star op- 
erators is equivalent to the product of all bulk star operators). 
For example, Gj^G^ is generated by the bulk star operators 
alone in bipartitions 2 and 3 in Fig.|2] while the product of all 
boundary star operators along one of the two boundaries must 
be included to generate G^G^ for bipartitions 1 and 4. 

Let us define B*" = {^''}, '6' for 'bulk', to be the set of Ising 
spin configurations on the sites of the square lattice where all 
9\ corresponding to boundary sites s are fixed to equal +1. 
Let us also define ~ {(^^}^ '5' for 'boundary', to be the 
set of Ising configurations where 9^ = +\ for all bulk star 



operators, 0f = +1 for all boundary star operators belong- 
ing to one chosen boundary, and 6'f = ±1 for the remaining 
boundary star operators, so long as all 9^ spins belonging to 
the same boundary have the same sign. Notice that {1} 
for bipartitions 2 and 3 in Fig. |2] where 1 is the configuration 
with all the spins 9^ — +1. One can finally show that there is 
a one-to-one correspondence between the elements of Gj^G^ 
and the Ising configurations in {^''0^ M 0^ e e\ 0* G 6*'}, 
where 0^6^ represents the site-by-site product of the two con- 
figurations (i.e., {9^9^)^ ^ 9^91). Therefore, 



E 



(42) 



and in particular. 



9l9l9M0A9)di0l, = 

{s,s') 

s , s bulk 

= 019^9 A9)0l- (43) 

{s,s') 

s' boundary 
s bulk 

+ 0l9M0A9)0i (44) 

s^s' different boundaries 

+ E 9i9,{g)9A9)ei (45) 

s,s' same boundary 

+ Y (^sigWs'ig), (46) 

{s,s') 

where we used the fact that if s is in the bulk then 6'f — +1, 
if s belongs to a boundary then 9^ ~ +1, and if both s and s' 
belong to the same boundary then 0f 0f/ = +1. 

Let us focus on the bipartitions of interest to compute the 
topological entropy ( |23] l. First of all, in the limit r,R ^ oo 
there are no nearest-neighboring stars s and s' belonging to 
two different boundaries. Therefore, the term i45[ vanishes 
identically. For bipartitions 2 and 3, 9* = {1} and 



The rh.s. of the above equation can be interpreted as the partition function of an Ising model with nearest-neighbor 
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interactions, where only the bulk degrees of freedom are 
allowed to flip starting from a given configuration 8{g). 
Clearly such partition function is invariant upon changing the 
initial configuration as long as the new one is in the same 
ergodic sector. 



For example, one can equivalently choose 

0(9)-- 

and the expression above simplifies to 



1 if s belongs to the bulk 
(g) if s belongs to the boundary, 



(48) 



Here ^23(5) represents the partition function of an Ising 
model with nearest-neighbor interaction of reduced strength 
J/T = (3, and with fixed spins along the boundary of bipar- 
titions 2 and 3, respectively. The values of the spins at the 
boundary are determined by g. 



For bipartitions 1 and 4, 8* = {l,f}, where the configura- 
tion f has all the spins equal to + 1 except for those belonging 
to the chosen boundary, say boundary 2 in Fig. |2l which are 
equal to —1. In this case 



^ .j-^s,s' sameboundiiry ^ . , , , ,^s,s'bulk^ ^ ^ s ' boundiiry 1 , s bulk a / \ 

heGj^ e''ee'' 
keOg 



= zlM + ztr {g)- (50) 

I 



Here Zi^{g) are the analog of Z 23(g) for bipartitions 1 

and 4, respectively, while '^'""^"^'^(g) differ from the for- 
mer by the fact that all the (fixed) spins belonging to bound- 
ary 2 in bipartitions 1 and 4 respectively have been flipped. 
In other words, Zf ^[g) represents the partition function of 
an Ising model with nearest-neighbor interaction of reduced 
strength J/T = (3, and with fixed spins along the bound- 



ary of bipartitions 1 and 4, respectively. The partition func- 
tions zf 4*'^"'''(g) differ in that the spins along one of the two 
boundaries have been flipped with respect to their values in 
Zf^{g). Again, the values of the spins at the boundary are 
determined by g. 

In this notation, the topological entropy of the system can 
be written as 



S, 



topo 



= lim 



9, twisted 



(g) Zf(g) + Zf"(.g) 



geo 



ziig)zi{g) 



(51) 



where the sum over g acts as a weighed average of the loga- 
rithmic term over all possible values of the spins at the bound- 
ary. Notice that in Eq. ( ISTT i the partitions with two boundaries, 
and hence with non-trivial topology, are those that appear 



with two contributions (bipartitions 1 and 4), corresponding 
to some relative boundary conditions (BCs) and their twisted 
counterparts. These contributions, as we show below in de- 
tail, are responsible for the non-vanishing topological entropy. 
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In the topological phase, the two partition functions for the 
twisted and untwisted BCs contribute equally, and in the non- 
topological phase, one partition function is exponentially sup- 
pressed when compared to the other, in the thermodynamic 
limit. Therefore, there is an extra entropy contribution in one 
of the phases depending on whether the boundaries of topo- 
logically non-trivial bipartitions are twisted or not relative to 
one another. 

From Eq. ( BTI ). the behavior of the topological entropy can 
be qualitatively argued as follows. Deep in the disordered 
phase, where the correlations are short ranged, the choice of 
boundary conditions is likely to affect the partition function of 
the system only with exponentially small corrections. Thus, 



we can expect to have Zf{g)Z^{g) ~ Z 



d. twisted 



(5)Zf (5) and 5, 



topo 



2. On the other hand, 
deep in the (ferromagnetically) ordered phase the partition 
function of a system with twisted boundary conditions is ex- 
ponentially suppressed with respect to the one without the 
twist. Thus, Zf{g) » Zf''™'^'^'^(g), Z^ig) » 
while Zf{g)Zi{g) ~ Zl{g)Z§{g) still holds. This leads to 

— n 

'-'topo ~ 

In the following two sections we will show with rigorous 
arguments that the behavior of the topological entropy across 
the transition is strongly first order, with a sudden jump from 

'5'topo = 2 to S'topo = 0. 



1. The disordered phase f/3 < (3^) 

In the limit of small f3, namely above the ordering transi- 
tion, one can compute S^^^^ via the high-temperature expan- 
sion of the Ising model with fixed spins at the boundary. 

Let us rewrite. 



^l(ff) = E n (cosh /3 + sinh/? 0,^0,, 

Qb^Qb bonds? 



2^/ 



2-SN 



Q seSg 



(52) 



where Sj , s'^ are the sites at the ends of bond i, N is the total 
number of bonds on the lattice, i5 is the length of the boundary 
in number of 9 spins, c = cosh (3, t — tanh f3. The sum over 
Q runs over all possible graphs on the bonds of the square 
lattice, composed entirely of closed loops and open strings 
connecting two boundary spins. The product HsGfi, ^sls) 
compasses all the boundary spins that appear as end points 
(the set £g) of open strings in Q. Finally, £{Q) is the total 
length of the closed loops and open strings in Q. 

Analogously for Zf (5). The case of Zf (5) and Zf (g) dif- 
fers from Z^ig) and Z^{g) in that there are now two types 
of open strings: those going from one boundary to itself and 
those connecting the two boundaries (see Fig.|2]l. 

Next, let us compare the product ((?)Zf (g) with the 
product Zf[g)Zf{g). Notice that bipartitions 2 and 3 have 
precisely the same total combined boundary as bipartitions 1 
and 4. In order for a graph to appear in one of the two products 



and not in the other, it needs to comprise loops or strings that 
are able to tell the difference between to two possible origins 
(2 H- 3 vs. 1 H- 4) of the total combined boundary. 

Examples of such open strings are shown in Fig. [3] One 





FIG. 3: (Color online) - Examples of open strings that appear in 
the expansion of Zf ((7)^3 {g) but are not present in the expansion of 
Zi [g)Z'l[g) (Top), and vice versa (Bottom). Top panel: The thick 
yellow and blue lines correspond to the boundaries in bipartitions 2 
and 3 respectively. Boundaries belonging to both are shown in a thick 
dashed yellow-blue pattern. The strings in question are symbolically 
represented by thin black lines. Bottom panel: Same color coding, 
with yellow corresponding to bipartition 4 and blue corresponding 
to biparition 1. Notice that strings appearing in one expansion and 
not in the other must connect boundaries of the same solid color, and 
therefore cannot be shorter than R~ 2r. 

can show that these telltale strings, and the analogous closed 
loops, cannot be arbitrarily short, and their length is bounded 
from below by i? — 2r. As a consequence, the correspond- 
ing graphs are exponentially suppressed at least as t^^"^"^ , and 
in the limit r, R ^ 00 with R — 2r 00, implicit in the 
definition of the topological entropy, one obtains 



^f(g)^f(.9) 
zi{g)zi{g) 



(53) 
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Similar considerations apply when comparing the product 
Ziig)Zi{g) with products of the kind zf'"^''""^{g)Zi{g). In 
this case, the boundaries involved are exaclty the same, and 
the relevant telltale elements of the graph are open strings con- 
necting one of the two components of the boundary with the 
other. Such strings are in fact the only elements that are sensi- 
tive to the twisted boundary conditions. Clearly the length of 
these strings is bounded from below by i? — 2r, and 



zf''-"=-'(g)Zf(g) 



(54) 



exponentially fast, at least as t^^"^"^ , with [R — 2r) — > oo. 

Of course, our reasoning is correct up to the point where the 
high temperature expansion breaks down, and entropic contri- 
butions balance the exponential suppression. Said differently, 
this is the case when the correlation length in the Ising model 
goes to infinity, and the large r, R limit does not guarantee that 
the ratios of products of partition functions above tend to one. 

Given Eqs. ( |53I54| |, we can finally use Eq. JSTT l to obtain the 
topological entropy of the system throughout the disordered 
phase P<I3^~ 0.4406868, 



see 
log2 4 = 



2. 



2. The Landau-Ginzburg ordered phase (f3 > (3^) 

What happens below this transition? Rather than attempt- 
ing a low-temperature expansion, it is convenient to use the 
duality relations derived by A. Bugrij and V. Shadura in 
Ref. 16 for the inhomogeneous, finite-size Ising model. In 
particular, they obtained the duality relations for a square lat- 
tice Ising model wrapped around a cylinder of finite length, 
with fixed, free and mixed boundary conditions. Follow- 
ing the usual convention, let us label /3 the coupling con- 
stant of the dual Ising model (defined on the plaquettes of 
the original lattice), which is related to (3 by the duality re- 
lation sinh/3 sinh/3 = 1. Let us also indicate with Z{6,6') 
and Z{S,6') the partition functions of the system on the fi- 
nite cylinder and its dual, with 6, S' specifying the boundary 
conditions, namely S, S' — O, x for free and fixed boundary 
spins, respectively. With this notation in mind, the results by 
Bugrij and Shadura - to the purpose of the present paper - can 
be summarized hy^ 

ZiO,0) - /C[Z(x,x) + Z'"'^"=''(x,x)] (55) 
Z{x,o) = /CZ(o, x) (56) 
Z(0, x) = /CZ(x,0), (57) 

where x , x ) differs from Z( x , x ) by the fact that the 

fixed boundary spins at one end of the cylinder have been 
flipped. Notice that the proportionality coefficient K, is the 
same in all the equations, and that Z{0, x ) = Z(x , O) and 

Z(x,o) = Z{0, x). 



Let us then consider in Eq. ( |49] ). Thanks to the nearest- 
neighbor character of the interaction between 6 spins, subsys- 
tem A interacts only with itself and with the boundary Jj, and 
so does subsystem B. Thus, one can factorize the two subsys- 
tems and obtain (cfr. Eq.(|49]l) 



E 




^ s.a' talk g,^,^ ^ b„„„d.„,. bult ^ 



(58) 



Similar arguments apply to bipartition 3, 



E 



heG^ 
keG, 



(59) 

and with a few more steps, to bipartitions 1 and 4 as well (cfr. 
Eq.®), 



E 

h€G^^ 
keG, p 



JJT.i<^Uh9k) ^ 



X e 



ZL{9) + Zfr"'\9)\ Zts^ [g) (5), 

(60) 



and 

E e 

kdG^ 



X 6 \ ■ ' 



y^Zl. (5) Zl.{g) ZU9) + 



(61) 
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where IB^ and IB2 refer to the two connected components 
of subsystem B in bipartition 1, i.e., the component inside 
boundary 1 and the component outside boundary 2, and anal- 
ogously for 4:A-^ and 4^2. 

In order to apply Eqs. ( 1551157b to the present case, some fur- 
ther considerations on the bipartitions in Fig. |2] are needed. 
Recall that, although S^^^^ is indeed a quantity of order one, 
we expressed it in Eq. (pTT i in terms of a ratio of extensive 
partition functions Zf{g). Thus, any sub-extensive correction 
to these partition functions (i.e., 0{2^ ), with a < 1, N be- 
ing the number of degrees of freedom in the system) will only 
amount to an exponentially small correction to S^^^^, that van- 
ishes in the thermodynamic limit. In this context, the partition 
function Zfj^lg) (see Fig.|4]l is 'equivalent', in the thermody- 
namic limit, to the partition function of an Ising model on an 
infinite cylinder with fixed boundaries at the edges (bound- 
ary 1 and 2, respectively). Similarly, the partition function 
Zfg [g) can be regarded as that of an Ising model on an in- 
finite cylinder with fixed boundary conditions on one edge 
(boundary 1 ) and open boundary conditions on a suitably in- 
troduced boundary 7]^. Finally, the same approach can be 
used for Zfg (g), with fixed boundary conditions on one edge 
(boundary 2) and open boundary conditions on another suit- 
ably introduced boundary 74. Qualitatively, this is illustrated 
in Fig. 2^, where the spins on boundary 1 and boundary 2 are 
fixed and those belonging to and 74 are free. Notice that the 





(b) 











nda: 


nda: 




3 


3 












XI 







(d) 



mapping onto infinite cylinders requires the distance between 
any of the 7^ boundaries introduced in Fig. |4] and any of the 
original boundaries in Fig.|2]to diverge with system size. Sim- 
ilar arguments apply to bipartitions 2, 3 and 4. This leads to 
a correspondence between our factorized partition functions 
and those used in Refill namely 



^i\(5),^i%(5),^2a(5), 
Zls{g),ZlAg),Zls{g). 
^4^(5), ^4% (5) 
zUg).ZtB{9) 

Zd, twisted / \ ryd, twisted / \ 



^(0,x) 
Z{x, x) 



\x,x). 



The results in Ref. fl^ can then be applied to our systems 
and lead to the following equations: 



rvd, twisted / \ 
ZlA id) 



yd 



(x,o) a Zfs^(g) 
,(x,o) « Zfeig) 



Z^a{x,o) (X 



^2B 



2A 

(X,0) CX ZIb 



(5) 
(5) 



Zf^(x,0) (X Zl^ig) 
Z|b(x,0) « z^sig) 



ZU0,0) (X ^fB(5) + ^4'B"(3) 



^4\(x,0) CX 
^fA,(x,0) CX 



(62a) 
(62b) 
(62c) 



(62d) 
(62e) 



(62f) 
(62g) 



(62h) 
(62i) 
(62j) 



For convenience of notation, let us define the dual partition 
functions for the whole system in the different bipartitions 



Zf = Zf^(0,0)Zi\(x,0)Zi%(x,0) (63) 

Z| = ^2a(x,0)Z|s(x,0) (64) 

^ Zf^(x,0)Zfs(x,0) (65) 

Zi EE Zt(x,0)z!.{x,0)Zis{0,0) (66) 



FIG. 4: (Color online) - Illustration of the four bipartitions used to 
compute the topological entropy in Ref. 10, with a possible choice 
for the additional boundaries needed to map each partition onto an 
infiiute cylinder in the thermodynamic limit. 



Finally, we have all the ingredients to evaluate the topolog- 
ical entropy for (3 > (3^. Let us first rewrite Eq. (ISTT l using 
Eqs. (ESMB instead of Eqs. (|49]|50li: 
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^ r,R^oo I Z 



(g) Zf„ (5)Zf. (g) Zf5(.g) + Z«''(g) 



^lA(5)^lB(5)^3A(ff)^li3(5) 

s,s'G(5]^+(54 s,s G(52+53 



(67) 



where the last two terms inside the square brackets come from the exponential factors in Eqs. ( I58H61I ). Using the duality 
relations ( |62] l. we can identify 



Zf^(5) + Z,T (5) Zf^(g)Zf^(g)Z!.ig)Z!^(g) Z^g) + Z: 



^d. twisted 



(5) 



Z^^A{9)ZlB{9)ZlA{9)ZtB{9) 

Zf^(0, 0)Zf„ (X , 0)Zf„ (X , 0)Zf^, (X , 0)Zf , (X , 0)^|b(0, O) 



7(3 yd ' 



ZlJx, 0)Z9^{x, 0)Zf^(x , 0)Zis{x, O) 



For (3 > /?^, the dual Ising models are in the disordered 
phase and one can perform a high-temperature expansion to 
calculate the ratio (Zf Zf ). Using the same loop 

description as for the original system, with t replaced by t = 
tanh(/3), one can show that (Zf Zf Z|) = 1 in the 
thermodynamic limit. 

The remaining terms in Eq. ( |67] | can be dealt with more 
promptly by reverting back to the original a spin degrees of 
freedom, 

s,s' €=(5^ +(5^ s,s'g52+S^ 

E OsigWAg) - E (^sigWAg) = 

{s,s') {s,s') 



= E E '^^^9): 



where i labels the bonds of the square lattice, and i £ 6 means 
that the bond i connects two sites s and s' belonging to S. This 
contribution can be shown to vanish identically since the set 
of boundary a spins in bipartitions 1 and 4 is identical to the 
set of boundary spins in bipartitions 2 and 3 (see Fig. (|2|i). 

In the end we find that 



5*, 



topo 



= lim 



sec 



ycf yd 
'^1 ^3 



= 



identically in the ordered phase P > (3^ 



B. Beyond 1-body potentials 

As we akeady mentioned, the calculations carried out in 
Sec. 1111 Al for the specific model presented in this paper can 
be straightforwardly extended to the case of any factorizable 
wavefunction. All one needs to do is identify a proper set of 
local generators (i.e., acting on the a spins contained within a 
disc of finite radius) for the group G, and the equivalent of the 
collective boundary flip operators. The rest of the derivation 
follows essentially unchanged, in the limit r,R^ 00. 

What happens if we attempt to generalize our approach fur- 
ther and we consider non-factorizable wavefunctions? For 
simplicity, take once again the Kitaev-like GS wavefunc- 
tion in Eq. ( fTSl l. but replace the argument of the exponential 
/? J2i '^iig^ Q;)/2 with some generic function ~(3Eg/2. As we 
can see immediately from Eq. ( l24b . our approach to compute 
the topological entropy can no longer be used from the very 
first stage. On the other hand, it is tempting to conjecture that, 
so long as Eg is short ranged (i.e., it can be written as the sum 
of terms involving cr spins within a disc of finite radius on the 
lattice), the error that one makes by neglecting the terms in- 
volving spins across the boundary of a bipartition (Eg) does 
not give topological contributions to S'topo- Under this assump- 
tion, one can then set E^ = and use the approximate equal- 
ity Eg ~ E^^ + Eg^ to re-establish the factorability needed 
to carry on with the calculations. The result obtained for S^^p^ 
in Eq. ( BTI ). which employed this approximation, nonetheless 
shows no explicit dependence on it in the final expression, and 
one could then reinstate the full Eg at that stage. If the conjec- 
ture above is correct, the formula in Eq. f5i\ gives the exact 
topological entropy for a generic GS wavefunction that satis- 
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lies (i) the positive amplitude condition, (ii) the group condi- 
tion for G,i2. and (iii) the locality (i.e., short ranged) condition 
on Eg. 

A simple example where this conjecture can be applied 
rather straightforwardly is the case where Eqs. (fl4b and (fTSl l 
are replaced by 

H = i?K.,aev + Ai^e-'^^.^=^(-)..««*-*^ (68) 

s 

\GS) = ^^e'^^(->"'(s)"^(s)/'<?|0). (69) 

(For a discussion of the general construction scheme of such 
type of Hamiltonians, see Ref. 13,) Here \ GS) is the GS wave- 
function of H in the topological sector where Yiiep '^f ^ 
The notation (ij) ,j^s stands for j nearest-neighbor of i but 
not adjacent to the same vertex s. Without loss of generality, 
we consider the range /3 S (0, oo), where the new term in the 
Hamiltonian favors ferromagnetic order in the a spins. 

Let us then introduce the same description in terms of the 
6 spins, as in the previous section. Given that the product of 
two nearest-neighboring a spins translates into the product of 
two nejcf-nearest-neighboring 9 spins, we obtain 



no local order parameter exists that resolves the topological 
phase. 

What is the fate of the topological entropy across this 
Landau-Ginzburg phase transition? According to the conjec- 
ture above, we can directly substitute ^^<j\ '^i^j 
into Eq. ( ISTI i and compute S^^^^. This amounts to replac- 
ing the boundary Ising partition functions in the argument of 
the logarithm with the partition functions of two decoupled 
boundary Ising models. As a result, all calculations carried 
out in the previous section remain essentially unchanged and 
one arrives to the identical result that S^^p^ — 2 throughout the 
high-temperature phase, and vanishes otherwise. 

This scenario is in agreement with previous results on a 
quantum version of the eight-vertex model by Ardonne et 
a/.,'** whose GS is a generalization of the one in our model. 
Consider indeed the wavefunction in Eq. ( |69] l. Given the na- 
ture of the group G, the four spins belonging to any plaquette 
of the square lattice can assume only eight distinct configura- 
tions (niep'^'i — +1)' illustrated in Fig.|5] Such configura- 
tions map naturally onto the vertices of an eight-vertex model 
upon replacing each positive spin with an arrow along the cor- 
responding bond of the dual lattice, pointing, say, from sublat- 
tice A to sublattice B, and vice versa for the negative spins (as 



g&G 



{ij)es ' ' ' 



the 



see 



(notice the additional factor of 2 in the exponent due to the 
fact that the same product 9^6^, corresponds to two distinct 
products crfaj), and 



\GS) = 



see 



Z 



-9i0)\0), 



(70) 



shown in Fig.|5]). Given that 

amplitudes in the GS wavefunction i69[ factorize into prod 
ucts of vertex fugacities 'i^p{f3J2{ij)£s '^f'^j Z^)- nota- 
tion of Fig.|5] the vertex fugacities assume the values a — b — 
1, c = e^^'^ and d = e^^ . The GS spatial properties of our 
model are therefore captured by a classical eight- vertex model 
with the appropriate fugacities'** (but see Ref. 13 for a general 
discussion of such quantum-to-classical correspondence), and 
one can then use Baxter's exact solution'^ to obtain the phase 
diagram as well as the scaling exponents at the critical point. 
All this is discussed in detail in Ref. the model under- 
goes a second-order, Z2-symmetry-breaking phase transition 



where Z = J^eee . The latter is the parti- when (P = + 2 (i.e., ^ (1/4) ln(V2 + 1)), separating 



tion function of a square-lattice Ising model with sole next- 
nearest-neighbor interactions, which factorizes into the prod- 
uct of two decoupled Ising models with nearest-neighbor in- 
teractions (namely corresponding to the 9 spins on each of the 
two sublattices). In this case, all equal-time correlators in the 
GS of the quantum system can be written in terms of classi- 
cal correlators of two decoupled Ising models. As before, we 
expect the system to undergo a phase transition when the two 
Ising models become critical at = (1/4) ln(V2 + 1) ~ 
0.2203434. However, contrarily to the previous case, the new 
model undergoes a spontaneous symmetry-breaking (Landau- 
Ginzburg) phase transition! This is best seen by mapping the 
system onto a quantum eight-vertex model, as discussed be- 
low. The local order parameter that captures the transition 
is the magnetization of the a spins, whose expectation value 
can be written as an ensemble average of the product of two 
neighboring 9 spins, i.e., belonging to two decoupled Ising 
models. Clearly such average vanishes identically in the high- 
temperature phase (3 < j3^, while it becomes finite in the 
ordered phase (3 > (3^. Notice that this local order param- 
eter that acquires an expectation value does so in the non- 
topologically ordered phase, as expected from the fact that 



a topologically ordered liquid phase from a Landau-Ginzburg 
ordered phase. The local order parameter across the transition 
is indeed the magnetization in the original a spins. 



IV. CONCLUSIONS 

In this paper we studied a topological quantum phase tran- 
sition in a microscopic model that can be examined analyt- 
ically. For this system, an extension of the toric code, the 
ground state wavefunction can be written exactly as a func- 
tion of the parameter (3 that drives the system across the quan- 
tum phase transition. We computed the topological entropy 
for this system as a function of (3, and showed that it remains 
at a constant non-zero value throughout the topologically or- 
dered phase (/3 < ~ 0.4406868). Immediately after the 
quantum phase transition at f3c, the topological entropy drops 
to zero and remains so in the non-topologically ordered phase 

The GS wavefunction of our quantum system has positive 
amplitudes in the basis of choice. This property allows us 
to relate many quantities that are relevant in characterizing 
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FIG. 5: Illustration of the eight allowed spin configurations around a plaquette of the square lattice, in any basis state g\0}, g £ G. These can 
in turn be mapped onto the configurations of an eight vertex model by replacing each positive spin with an arrow along the corresponding bond 
of the dual lattice (dashed lines), pointing, say, from sublattice A to sublattice B (vice versa for the negative spins). Clearly, the corresponding 
eight vertex configurations differ depending on the location of the plaquette in the dual lattice: sublattice A (middle), or sublattice B (bottom). 
The letters a, b, c, d correspond to the usual labeling of the vertex fugacities in the eight vertex model. 



the (2+1 D) quantum system to those of a simple (2D, not 
3D) classical Ising model at an inverse (classical) temperature 
equal to the value of the coupling constant /? that drives the 
quantum system through the T = phase transition. For ex- 
ample, the magnetization of the quantum system equals the 
energy EisingiP) of the classical Ising model. While the mag- 
netization is continuous and non-vanishing across the quan- 
tum phase transition (much as the energy i^ising (l^) is across 
the classical Ising transition), its derivative with respect to /3 
diverges logarithmically at f3c (much as the Ising model heat 
capacity Cising diverges logarithmically at /3c). 

Despite the relation to the 2D classical Ising model, the 
quantum phase transition does not have a local order parame- 
ter that vanishes on one side and not on the other. Of course 
one expects that no local order parameter can characterize the 
topological phase, but in this particular example, there is no 
order parameter that characterizes the non-topologically or- 
dered phase either. One can indeed identify from the map- 
ping to the Ising model a parameter that orders in the non- 
topological phase; however, this variable is non-local in the 
physical spin variables used to define the local Hamiltonian. 
Specifically, the order parameter is, in the language used in 
this paper, the expectation value of the Og variables defined 
on the sites of the square lattice, such that cr(ss/) = SgO'^, for 
nearest neighboring sites s, s' . While cr^s^') is obviously local 
in terms of the 9s;0'^, the inversion needed to write the 6''s in 
terms of the cr's is non-local. Hence, {Os) may detect the tran- 
sition into the non-topological phase, but as it is non-local it is 
not an order parameter in the usual sense. That there is no or- 
der parameter for the non-topologically ordered phase is not 
generic (see the example in Sec. IIIIBI ), as perhaps the most 
obvious exit from a topological phase is by escaping into a 
locally ordered phase due to spontaneous symmetry breaking. 
Hence, the main example studied in this paper is particularly 



interesting in that one has no local order parameter in either 
phases. 

Recently, P. Zanardi et al. proposed a new approach to 
study quantum phase transitions through the behavior of a 
metric tensor in parameter space, derived from a fidelity-based 
notion of distance between states.^" In this approach, one does 
not require any a priori knowledge of an order parameter to 
detect a phase transition (one may argue that the fidelity en- 
compasses correlations of both local and non-local operators). 
Thus, these ideas may be particularly useful to detect topolog- 
ical quantum phase transitions (see Ref. 9). 

We end with a speculative note: topological quantum phase 
transitions should share the feature that, even in the case when 
there is no local order parameter in either the topological or 
the non-topological phase, there should be, generically, de- 
tectable singularities in high enough derivatives of local ob- 
servables with respect to the coupling that takes the system 
across the transition (as in the case study presented here)i^ 
After all, what would be a phase transition without singulari- 
ties in any physical observable? 

Upon completion of this work, we became aware of simi- 
lar work being pursued from a gauge theoretical perspective 
by S. Papanikolaou, K. Raman and E. Fradkin on the quan- 
tum eight-vertex model,^^ to which we refer the reader for a 
complementary approach. 
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